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1.  INTRODUCTION 


Lidar  backscatter  from  cirrus  cloud  particles  at  multiple  wavelengths  was  measured  on 
several  days  at  Hanscom  AFB,  Massachusetts  in  August  2001 .  The  range-resolved 
measurements  were  made  for  several  hours  each  day  at  altitudes  from  9  km  to  16  km  over 
a  temperature  range  of -30  C  to  -  70C.  The  wavelengths  of  the  lidar  measurements  were 
0.532  pm,  1.064  pm,  and  10.6  pm.  Lidar  backscatter  provides  useful  information  on  the 
vertical  structure  of  cirrus  and  its  temporal  variation.  Backscatter  measurements  are 
analyzed  using  scattering  theory  calculations  to  better  understand  their  behavior  with 
respect  to  wavelength  and  characteristic  cirrus  particle  sizes.  The  backscatter  also 
yields  information  on  microphysical  and  optical  properties  of  cirrus  such  as  number 
density  (ice  water  content),  extinction,  and  extinction-to-backscatter  ratio,  which  can  be 
useful  in  laser  propagation  and  atmospheric  science  applications.  This  report  describes 
these  measurements,  relevant  calculations,  and  their  comparison  with  previously- 
published  results. 


2.  BACKSCATTER  FOR  SINGLE,  INDEPENDENT  SCATTERERS 


Backscatter  is  the  efficiency  of  an  object  for  scattering  electromagnetic  radiation  back 
toward  the  source  at  the  scattering  angle  8  =n.  Figure  1  (after  van  de  Hulst,  1957) 
depicts  light  from  below,  incident  on  a  scattering  particle  and  an  associated  scattering 
angle.  Backscatter  is  quantitatively  represented  as  the  scattering  cross  section  for 
scattering  angle  6  =  n  divided  by  the  object’s  geometrical  cross-section.  For  a 
homogeneous  sphere,  the  backscatter  efficiency  Qbscl  can  be  expressed  in  terms  of  the 


scattering  amplitude  function  5,  as 


Qbsct 


(1) 


1 


lA 


it  A  A  A  a 


Figure  1 .  Incident  light  from  below  on  a  scattering  particle  or  object.  6 
represents  the  scattering  angle,  (after  van  de  Hulst,  1957) 


(Case  I)  (Case  II) 


Figure  2.  Incident  light  from  below  on  a  long  cylinder.  6  represents  the 
scattering  angle  and  a  is  the  angle  of  obliquity,  (after  van  de  Hulst,  1957) 


where  St  is  the  scattering  amplitude  function  evaluated  at  0  =  n  and  x  =  kr  =  ——  is  the 

A 

particle  size  parameter  (van  de  Hulst,  1957).  Here,  X,  is  the  wavelength  of  incident 
radiation  and  r  is  the  sphere  radius.  The  scattering  amplitude  5,  is  a  function  of  the 

incident  energy  wavelength  A ,  particle  size  or  radius  r,  and  particle  refractive  index  m, 
which  is  a  complex  number.  For  a  collection  of  spheres  of  different  sizes,  i.e.  a 
polydispersion,  the  volume  backscatter  coefficient  is  written  as 

Pn  _  Qhsct n(r) nr1  dr  (m'1  sr'1)  _  (2) 

An 

■  •  -3 

where  n  (r)  is  the  particle  number  density  size  distribution  (the  number  of  particles  m 

pm'1).  The  variables  shown  in  Eqs.  (1)  and  (2)  will  be  described  further  in  this  report 
in  addition  to  the  scattering  geometry  associated  with  light  incident  from  below  on  a  long 
cylinder,  as  depicted  in  Figure  2. 


3.  LIDAR  CHARACTERISTICS  AND  MEASURED  LIDAR  BACKSCATTER 


Lidar  provides  measurements  of  range-gated  or  range-resolved  backscatter  from 
scattering  particles  at  optical  wavelengths.  For  the  measurements  described  in  this  report, 
two  lidar  systems  were  used  to  provide  the  range-resolved  backscatter.  The  first  system 
employed  was  a  Nd:YAG  lidar  capable  of  backscatter  measurements  at  three 
wavelengths  (0.355  pm,  0.532  pm,  and  1.064  pm).  However,  only  two  wavelengths, 
0.532  pm  and  1.064  pm,  were  operating  during  the  described  measurements.  A  second 
lidar,  a  C02  heterodyne  detection  lidar,  was  used  and  made  range-resolved  backscatter 
measurements  at  10.591  pm.  The  two  lidar  systems  were  physically  adjacent  to  one 
another  during  these  measurements.  Table  1  lists  the  operating  characteristics  of  each 
lidar  system. 
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Table  1.  Lidar  Characteristics 


Lidar  Nd:YAG  C02 


Wavelengths  (pm) 

0.532,  1.064 

10.591 

Telescope  aperture  diam. 

30  cm 

30  cm 

Range  gate 

120  m 

120  m 

Pulse  rep.  frequency 

30  Hz 

30  Hz 

Pulse  integration  period 

30  sec 

30  sec 

Pulse  Energy 

250  mj  (0.532  pm) 

400  mJ  (1.064  pm) 

60  mJ 

Nominal  meas.  range 

>  100  km 

60  km 

Full  hemispheric  field-of-view  scan  (both  lidars) 


While  these  lidars  can  take  measurements  at  any  azimuth  or  elevation  angle  above  0°, 
most  of  the  lidar  cirrus  backscatter  measurements  were  usually  taken  at  near-zenith 
elevation  angles.  The  exceptions  were  for  most  contrail  measurements  and  for  some 
distant  cirrus  shields,  where  the  lidar  pointing  between  systems  was  generally  co-aligned. 
However,  these  off-vertical  axis  measurements  are  not  analyzed  in  this  report.  In 
addition  to  the  lidar  measurements,  radiosondes  were  launched  near  the  lidars  every  3 
hours  during  the  measurement  period.  The  radiosondes  measured  vertical  profiles  of 
temperature,  relative  humidity,  pressure,  and  horizontal  wind  velocity  from  the  surface  to 
20  km  above  sea  level. 


The  lidar  volume  backscatter  coefficient  is  calculated  from  the  C02  heterodyne  lidar 
measurements  by 


PAz)  =  kcar^~^  (nC'sr-')  (3) 

Lp 

where  kca;  is  the  absolute  calibration  coefficient  (s  tri1  sr ;),  Pr  is  the  range-  and 
extinction-compensated  power  (J  s'1),  and  Ep  is  the  transmitted  pulse  energy  (J).  For  the 
Nd:YAG  lidar  measurements,  the  measured  lidar  backscatter  is  determined  as 

Ar  (Z)  =  Pn Rayleigh  '  (clear)  ^  ^  ^ 
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where  PxRayleigh  is  the  Rayleigh  volume  backscatter  coefficient  for  air  molecules  and 
Pr  (clear)  is  the  range-compensated  power  in  the  clear  air  immediately  at  the  base  of  the 
cloud.  Rayleigh  backscatter  from  air  molecules  is  a  function  of  wavelength  and  air 
density  and  varies  as  A'4 .  At  cirrus  altitudes  (10  km),  typical  Rayleigh  backscatter 
values  are  ~  7  x  10'7  m'1  sr'1  at  0.532  pm  and  ~  6  x  10'8  m'1  sr'1  at  1.064  pm.  The 
measured  lidar  power  return  Pr  was  corrected  for  range  diminishment  using  the  1/r2 

formula  and  compensated  for  extinction.  The  extinction  compensation  was  accomplished 
using  a  vertical  profile  of  laser  line  molecular  transmittance  for  each  laser  wavelength 
computed  from  the  HITRAN  data  base  (Rothman  et  al,  1992)  and  radiosonde  input  data 
(pressure,  temperature,  and  relative  humidity)  nearest  in  time  to  the  measurements. 

Using  these,  a  calibrated  backscatter  profile  was  produced.  While  the  Nd:YAG  lidar  is 
self-calibrating  using  the  Rayleigh  backscatter  signal,  the  calibration  coefficient  kcai  re¬ 
quired  for  the  CO2  lidar  measurements  was  determined  from  time  integrated  target  par¬ 
ticle  board  backscatter  measurements  and  laboratory  analysis  of  the  particle  target  board 
samples.  For  both  the  Nd:YAG  and  C02  lidar,  the  returned  power  is  integrated  over  a  30 
sec  period  to  increase  the  sensitivity  of  the  signal.  For  the  C02  heterodyne  system, 
typical  measured  backscatter  magnitudes  from  cirrus  are  ~  10‘9  m'1  sr'1  while  typical 
cirrus  backscatter  magnitudes  from  the  NdrYAG  lidar  at  0.532  pm  are  ~  1 0 6  m 1  sr 1 . 

There  are  several  factors  that  can  affect  the  lidar  backscatter  measurement  and  its 
interpretation.  For  C02  lidar  backscatter,  the  accuracy  of  the  calibration  coefficient  is 
probably  only  within  +/-  20%  of  its  true  value.  Non-linear  effects  or  saturation  of  the 
power  return  signal  also  must  be  considered  for  both  lidars.  For  both  lidar  systems, 
specular  reflection  from  flat  cirrus  particles  or  plates,  which  can  produce  jumps  of  two  to 
three  orders  of  magnitude  in  the  measured  backscatter,  can  usually  be  avoided  by  tilting 
the  lidar  beam  slightly  off  vertical.  In  comparing  backscatter  magnitudes  at  different 
wavelengths,  changes  in  particle  composition  or  particle  morphology  within  the  cloud 
during  the  measurements  must  also  be  considered. 


4.  MEASUREMENT  LOCATION  AND  SYNOPTIC  WEATHER  CONDITIONS 


The  lidar  and  radiosonde  measurements  described  in  this  report  were  made  at  Hanscom 
Air  Force  Base  which  is  located  in  northeast  Massachusetts  at  latitude  42°  27'  N , 
longitude  71°  16'  W  at  70  m  above  sea  level.  Lidar  measurements  were  made  on  nine 
days  covering  the  period  7  August  2001  -  30  August  2001 .  The  measurement  site 
experienced  air  mass  types  of  different  origin  (maritime  tropical,  Canadian,  etc.)  during 
this  period  which  resulted  in  changes  in  the  tropopause  height  between  days”  Cirrus  was 
observed  during  the  nine  days  associated  with  deep  convection,  weak  waves  in  the 
synoptic-scale  westerly  flow,  and  surface  frontal  zones  and  jet  streams.  The  daily 
measurement  periods  generally  ran  from  1700  UTC  -  2400  UTC. 

Three  measurement  days  were  selected  for  study  -  9  August,  14  August,  and  16  August. 
These  were  days  where  extended  periods  of  relatively  uniform  cirrus  occurred  that  were 
not  associated  with  deep  convection.  On  these  days,  temperatures  at  the  cirrus 
measurement  altitudes  ranged  from  -30C  to  -70C  and  the  tropopause  height  ranged  from 
~  13  km  to  15.  5  km. 

On  9  August,  the  measurement  area  was  dominated  by  a  very  warm,  humid  air  mass  in 
the  troposphere  with  weak,  westerly  zonal  flow  at  upper  levels.  This  is  depicted  in 
Figure  3  which  is  a  contour  plot  of  the  geopotential  height  and  wind  speed  on  the  300  mb 
pressure  surface.  Extensive  cirrus  was  observed  from  1 300L  -  1 800L.  The  cirrus¬ 
generating  mechanism  was  likely  associated  with  a  weak  synoptic  scale  wave  in  the 
upper  tropospheric  westerly  flow.  With  a  tropopause  located  near  1 5.5  km,  the  cirrus 
was  measured  at  altitudes  from  12  to  15.5  km.  On  14  August,  a  Canadian  air  mass 
dominated  the  region  and  a  relatively  strong  polar  jet  stream  axis  oriented  SW  -  NE  was 
anchored  over  eastern  Massachusetts  (see  Figure  4)  with  peak  wind  speeds  exceeding  50 
m/s  from  9  km  to  13  km  altitude  as  measured  by  radiosonde.  A  distinct  cirrus  shield  was 
associated  with  the  jet  stream  and  the  shield’s  trailing  edge  lay  just  west  of  Hanscom 
AFB.  The  tropopause  was  lower  on  14  August  compared  to  9  August  with  the 
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Figure  3.  300  mb.  Analysis  10  August/OOZ  2001.  The  solid  lines  depict  geopotential 
height  contours  at  300  mb  (units  of  decameters).  The  dashed  lines  denote  wind  speed  in 
knots.  Hanscom  AFB  is  denoted  by  the  identifier  BED.  Concentric  circles  show 
location  of  synoptic  radiosonde  sites  used  to  compute  area-averaged  vertical  velocity. 
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Figured  300  mb.  Analysis  15August/00Z  2001  The  solid  lines  depict 
geopotential  height  contours  at  300  mb  (units  of  decameters).  The  dashed  lines 
depict  wind  speed  in  knots.  Hanscom  AFB  is  denoted  by  the  identifier  BED. 
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Figure  5.  300  mb.  Analysis  17  August/OOZ  2001  The  solid  lines  depict  geopotential 
height  contours  at  300  mb  (units  of  decameters).  The  dashed  lines  depict  wind  speed 
in  knots.  Hanscom  AFB  is  denoted  by  the  identifier  BED. 
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tropopause  altitude  being  13  km.  On  16  August,  the  polar  jet  stream  axis  had  moved 
well  east  of  Massachusetts  (see  Figure  5)  along  with  its  associated  cirrus  shield. 

However,  an  extensive  area  of  scattered  to  broken  cirrus  ahead  of  a  cold  front  in  the 
Upper  Midwest  advanced  into  New  England  and  persisted  during  most  of  the 
measurement  period.  The  area  remained  under  the  influence  of  a  modified  Canadian  air 
mass  and  the  tropopause  altitude  increased  slightly  to  1 3.5  km.  No  jet  stream  structure 
was  observed  and  weak  westerly  flow  dominated  through  the  troposphere  with  wind 
speeds  remaining  at  or  below  20  m/s. 

Three  daily  samples  of  time  vs.  altitude  profiles  of  uncalibrated  lidar  backscatter  imagery 
for  10.591  pm  wavelength  are  shown  in  Figures  6,  7,  and  8.  The  imagery  uses 
backscatter  profiles  which  are  produced  approximately  every  30  seconds.  The  light 
colored  areas  denote  backscatter  higher  than  the  detector  noise  floor.  The  areas  of 
enhanced  backscatter  below  4  km  are  associated  with  haze  or  aerosol  particles. 

Figures  9,  10,  and  1 1  are  plots  of  temperature  and  relative  humidity  as  measured  by  the 
Hanscom  AFB  radiosonde  launches.  Area-averaged  vertical  velocity  profiles  are  also 
plotted,  which  are  diagnosed  from  00Z  radiosonde  profiles  forming  a  triangle  -  Albany, 
New  York  (ALB),  Chatham,  Massachusetts  (CHH),  and  Gray,  Maine  (GYX)  -  whose 
sides  enclose  Hanscom  (BED)  (see  Figure  3).  Vertical  velocity  is  calculated  from  this 
triangle  using  the  kinematic  method  (O’Brien,  1970).  These  radiosonde  launches 
occurred  approximately  two  to  four  hours  after  the  cessation  of  daily  lidar  measurements. 

On  9  August,  the  tropopause  height  was  approximately  15  km,  the  base  of  the  cirrus 
backscatter  was  evident  at  ~  12.5  km  and  the  top  of  the  cirrus  backscatter  was 
approximately  1 5  km.  This  corresponds  to  a  temperature  range  of -50C  to  -70C  (see 
Figure  9).  The  tropopause  height  also  serves  to  cap  the  upper  altitude  of  the  cirrus.  On 
14  August,  the  tropopause  height  was  about  13  km  (see  Figure  10)  and  the  cirrus 
backscatter  extended  from  9  km  to  13  km  (Figure  7).  This  corresponded  to  a  temperature 
range  of  approximately  -30C  to  -60C.  Note  the  vertical  striations  of  cirrus  backscatter 
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Figure  6.  Time  vs.  altitude  Cross-section  of  Lidar  Backscatter  Return  for  10.6  pm 
for  9  August  2002.  Cirrus  clouds  are  observed  above  12  km.  Lighter  colored  areas 
below  7  km  are  associated  with  returns  from  haze  particles  and  water  droplet  clouds. 
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Figure  7.  Time  vs.  Altitude  Cross-section  of  Lidar  Backscatter  Return  for 
10.6  |im  for  14  August  2002.  Cirrus  clouds  are  observed  from  about  9  km 
to  1 1 .5  km.  Lighter  colored  areas  below  5  km  are  associated  with  returns 
from  haze  particles. 
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Figure  8.  Time  vs.  Altitude  Cross-section  of  Lidar  Backscatter  Return 
for  10.6  pm  for  16  August  2002.  Cirrus  clouds  are  observed  from  about 
9  km  to  13  km.  Lighter  colored  areas  below  7  km  are  associated  with 
returns  from  water  droplet  clouds  (cumulus  and  altostratus)  and  haze 
particles  below  3  km. 
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Figure  10.  Radiosonde  temperature  (C)  and  Relative  Humidity  (%)  Profile  14  Aug  19  GMT,  at 
Hanscom  AFB,  MA.  On  the  right  is  the  diagnosed  area-averaged  vertical  velocity  (cm  s'1)  for  1 5 
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Figure  11.  Radiosonde  Temperature  (C)  and  Relative  Humidity  (%)  Vertical  Profile:  16  August  19 
GMT  at  Hanscom  AFB,  MA.  On  the  right  is  the  diagnosed  area-averaged  vertical  velocity  (cm  s'1)  for 
17  Aug  00Z 


which  are  prevalent  on  this  day  but  absent  on  the  other  two  days.  On  1 6  August,  the 
tropopause  height  was  —  13.8  km  (Figure  11)  and  the  cirrus  altitudes  varied  from  9  km  to 
13.5  km  (Figure  8).  These  figures  illustrate  qualitatively  the  variation  in  cirrus  height 
and  structure  that  occurred  during  the  measurement  period. 

A  weak  qualitative  correlation  between  the  altitudes  of  maximum  relative  humidity, 
diagnosed  vertical  velocity,  and  altitudes  of  cirrus  backscatter  was  observed.  On  9 
August,  the  vertical  velocity  peaked  near  14  km,  in  the  cirrus  backscatter  region,  at  ~  0.7 
cm/s,  but  the  relative  humidity  was  measured  <  40%  in  the  upper  troposphere  region.  On 
14  August,  the  correlation  of  both  relative  humidity  and  vertical  velocity  with  cirrus 
backscatter  was  better.  The  vertical  velocity  peaked  at  >  10  cm/s  near  10  km,  and  a 
relatively  moist  layer  with  peak  relative  humidity  ~  70%,  was  observed  from  9  km  to  13 
km.  This  corresponds  closely  to  the  10.6  pm  cirrus  backscatter  vertical  extent.  On  16 
December,  the  relative  humidity  peaked  at  80  %  near  7  km  while  the  vertical  velocity 
showed  a  broad  peak  centered  at  8  km  with  a  maximum  of  4  cm/s.  These  peaks  were 
somewhat  below  the  cirrus  base  altitudes  of  ~  9  km  and  may  have  been  associated  with 
the  lower  altostratus  or  altocumulus  (see  Figure  8). 
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5.  CALCULATION  OF  BACKSCATTER  COEFFICIENT  FROM 
SCATTERING  THEORY 

Using  Eqs.  (1)  and  (2),  the  theoretical  volume  backscatter  coefficient  can  be  calculated 
for  a  polydispersion  of  spheres.  The  backscattering  amplitude  5,  (n)  is  calculated,  using 
the  expression  described  in  van  de  Hulst  (1957),  as 

-5,  (*)  =  S2  (*r)  =  X(n  +  y2)(-l)n  (a„  - K )  (5) 

n-\ 

where  an  and  bn  are  scattering  coefficients  for  spheres  and  are  functions  of  sphere  radius, 

electromagnetic  wavelength,  and  the  particle’s  complex  refractive  index.  The  scattering 
coefficients  are  calculated  using  Aden’s  (1951)  algorithm  employing  spherical  Bessel 
functions  and  Infeld’s  (1947)  recursion  formula  for  logarithmic  derivatives.  Upward 
recursion  is  used  for  spherical  Bessel  functions  with  real  arguments  while  downward 
recursion  is  used  for  those  functions  with  complex  arguments  (e.g.  described  by  Kattawar 
and  Plass,  1967).  The  wavelength-dependent  refractive  index  for  ice  is  determined  from 
the  tables  of  Warren  (1984).  Table  II  lists  the  complex  refractive  index  for  ice  for  the 
three  lidar  measurement  wavelengths.  A  description  of  the  calculation  procedures  for  the 
spherical  scattering  coefficients  an  and  bn  is  provided  in  Appendix  A. 

Table  2.  Refractive  index  for  ice  m  =  mr  —  i  mi  /  =  V-T 


Wavelength  ( /i/w  ) 

mr 

m, 

10.6 

1.1038 

0.1246 

1.064 

1.3004 

1.9  x  10'6 

0.532 

1.3117 

2.6  x  10’9 

Cirrus  clouds  are  typically  composed  of  non-spherical  particles  such  as  hollow  columns, 
bullet  rosettes,  dendrites,  hexagonal  plates,  capped  columns,  etc.  Using  ice  spheres  to 
approximate  hexagonal  ice  crystals  for  their  scattering  and  absorption  properties  can  be 
inadequate  and  frequently  misleading  (Takano  and  Liou,  1995).  In  an  attempt  to  address 
the  non-spherical  shape  of  cirrus  particles,  scattering  calculations  were  also  performed  for 
very  long  ice  cylinders  where  the  length  of  the  cylinders  is  significantly  larger  than  the 
diameter.  Here,  the  light  on  the  cylinders  is  assumed  to  be  of  perpendicular  incidence 
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from  below  and  the  cylinders  are  uniformly  oriented  in  the  horizontal  plane  such  that  the 
angle  of  obliquity  a  =  0  (see  Figure  2).  Following  van  de  Hulst  (1957),  the  backscatter 
efficiency  (per  unit  cylinder  length)  for  light  of  perpendicular  incidence  on  uniformly- 
oriented  circular  cylinders  can  be  expressed  as 

Case  I  (£  ||  axis):  Qbscr=—  ft  (*)f  =—*,(*)  ^  (6) 

TV  X  TV  X 

where  the  electric  field  E  is  parallel  to  the  cylinder  axis  and 

Case  2  (H  ||  axis):  Qbsct  =  —  \T2  (xf  =  —i2  (x)  (sc1)  (7) 

7VX  71 X 

where  the  magnetic  field  H  is  parallel  to  the  cylinder  axis.  The  scattering  angle  6  =  x , 
and  T}  and  T2  are  the  amplitudes  of  scattered  light  vibrating  parallel  and  perpendicular  to 

the  electric  field  vector,  respectively.  The  total  backscatter  efficiency  for  incident  light 
is  the  average  of  Eqs.  (6)  and  (7)  (see  van  de  Hulst,  1957)  or 

Qbsc =o.5f— i,  (*)+— i2  Mh  M)  (sr'‘)  (g) 

\nx  xx  J  xx 

Here, 

7U0)  =  6o  +  2j^cosn0  (9> 

«=1 

and 

T2(0)  =  aQ  +  2^jan  cosnO  (10) 

n- 1 

where  an  and  bn  are  the  cylindrical  scattering  coefficients  and  are  a  function  of  cylinder 
radius,  incident  wavelength,  and  the  scatterer’s  complex  refractive  index.  Their  solution 
is  expressed  in  terms  of  the  Bessel  functions  of  integer  order  (e.g.,  Kerker,  1969)  where 
forward  and  backward  recursion  techniques  can  be  employed  using  Infeld’s  (1947) 
logarithmic  recursion  relation.  A  full  description  of  calculation  procedures  for  the 
scattering  coefficients  for  a  long  cylinder  is  contained  in  Appendix  B. 

Analogous  to  Eq.  (2),  for  a  polydispersion  of  long,  finite  length  circular  ice  cylinders 
with  variable  radius  r ,  the  backscatter  coefficient,  per  unit  cylinder  length  /,  becomes 
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(nT2  sr'1) 


(ID 


In  the  evaluation  of  Eqs.  (2)  and  (1 1),  a  modified  gamma  particle  size  distribution 
(Intrieri  et  al,  1993)  is  used  to  represent  the  cirrus  cloud  polydispersion  and  is  of  the 
form: 

N  r  Y"'  ( 

=  X  —  exp  (# m'3  pm’1)  (12) 

r(wArJ  l  rj 

where  N0  is  the  number  density,  T  is  the  gamma  function,  p  is  an  integral  dispersion 
factor  and  rm  is  the  characteristic  radius  of  the  distribution.  The  characteristic  radius  rm 
of  the  size  distribution  is  related  to  the  equivalent  radius  re  by  re  =  (p  +  2)  rm .  For  the 
measurement  comparisons  with  theory,  p  =  2.  An  example  plot  of  the  modified  gamma 
size  distribution  n(r)  normalized  byN0  is  depicted  in  Figure  12.  Equation  (2)  is 

integrated  over  the  particle  radius  size  range  of  ~  0.01  pm  to  6  mm  for  spheres  and  Eq. 

(1 1)  is  integrated  over  the  cylinder  radius  size  range  of  0.01  pm  to  2  mm.  The  relation 
between  these  ranges  of  particle  size  radius  r  and  size  parameter  x  as  a  function  of 
wavelength  is  listed  in  Table  3  since  the  integrals  involving  scattering  efficiencies  are 
typically  evaluated  in  terms  of  x. 

Table  3.  Numerical  integration  limits  for  size  parameter  x  as  a  function  of 
wavelength  and  particle  type  x  =  Inrfk 

Wavelength  {pm)  Spheres  (xmin,xnax )  Long  Cylinders  (x^x^ ) 

10.6  0.006,3557  0.006,  1200 

1.064  0.059,35431  0.059,  11810 

0.532  0.118,70863  0.118,  23620 

In  the  backscatter  and  extinction  ratio  comparisons  to  be  discussed  subsequently,  the 
number  density  N0 ,  which  is  a  constant  in  Eq.  (12),  will  cancel  out.  The  number  density 

N0  can  also  be  estimated  from  the  expressions  defining  backscatter  and  extinction 
coefficients  for  a  polydispersion. 
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Modified  Gamma  Particle  Size  Distribution 


Figure  12.  Normalized  Modified  Gamma  Particle  Size  Distribution  Calculated  from 
Eq.  (12)  Used  to  Represent  the  Cirrus  Cloud.  Here  the  equivalent  radius  re  =  5  pm  and 
the  integral  dispersion  factor  p  =  2 
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6.  COMPARISON  OF  MEASURED  AND  CALCULATED  LIDAR  CIRRUS 
BACKSCATTER  RATIOS 


Intrieri  et  al  (1993)  compared  measured  and  theoretical  calculations  of  backscatter 
coefficient  ratios  for  lidar  and  radar  to  infer  the  effective  particle  size  of  cirrus.  A  similar 
approach  is  attempted  here  using  lidar  volume  backscatter  coefficient  ratios  for  the  three 
days  of  measurement:  9  August,  14  August,  and  16  August. 


The  backscatter  ratio  yn  for  two  lidar  wavelengths  given  simultaneous  measurements  at  a 
particular  altitude  z  and  time  t  is  defined  as: 


.Ml 

AM 


(13) 


where  A,  and  \  are  the  two  lidar  wavelengths. 


For  this  approach,  vertical  profiles  of  measured  lidar  backscatter  at  each  of  the  three 
wavelengths  are  interpolated  to  equal  altitude  levels  separated  by  100  m.  The 
interpolated  data  are  then  smoothed  using  a  Gaussian  probability  distribution  with  a 
standard  deviation  of  400  m.  Except  for  possible  missing  data,  backscatter 
measurements  from  the  Nd:YAG  lidar  at  1.064  pm  and  0.532  pm  are  considered 
simultaneous  in  time  and  altitude  location.  If  the  10.6  pm  (CO2)  and  0.532/1.064  pm 
(Nd:YAG)  lidar  wavelengths’  backscatter  measurements  are  within  16  seconds  of  each 
other  and  the  interpolated  height  levels  are  equivalent,  then  the  backscatter  ratio  yn  (Eq. 
(13))  is  evaluated  at  that  altitude  z  and  time  t. 

Measured  and  theoretical  backscatter  ratios  are  calculated  for  two  different  wavelength 
pairs:  10.591  pm/0.532  pm  and  1.064  pm/0.532  pm.  Lidar  backscatter  measurements  are 
used  from  three  days  in  the  calculations:  9  Aug  2001, 14  Aug  2001 ,  and  16  Aug  2001 . 
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Time  vs.  height  cross-sections  of  measured  logio  Pn  (m'1  sr ')  are  shown  in  Figures  13 

through  15  for  0.532  jam,  10.6  pm,  and  1.064  pm,  respectively.  These  may  be  useful  in 
certain  applications  by  providing  some  representative  or  characteristic  time  and  altitude 
variation  of  cirrus  backscatter  magnitudes  at  visible  and  infrared  wavelengths.  Figure  13 
displays  contours  of  the  logarithm  (base  10)  of  the  measured  volume  backscatter 
coefficient  at  0.532  pm  plotted  as  a  function  of  time  and  altitude  for  9  August  2001 

for  altitudes  above  10  km.  At  0.532  pm,  logio(^)  magnitudes  range  from  ~  -  4.5  to 
almost  -  6.0.  Figure  14  is  the  same  type  of  analysis  as  Figure  13  but  for  logio  P„  at  10.6 
pm  .  For  this  wavelength,  the  range  of  magnitudes  runs  from  ~  -6.5  to  -8.5.  For  1.06  pm 
(Figure  15)  on  9  August,  the  range  of  measured  log)0  {PK)  is  -5.0  to  -6.0. 

Figure  16  shows  contours  of  the  measured  backscatter  ratio  yn  for  10.6  pm  and  0.532 
pm  corresponding  to  the  cloud  backscatter  plots  of  Figures  13  and  14.  The  measured  yn 

ratios  range  from  about  0.003  to  0.01 .  These  can  be  compared  with  the  theoretical  ratio 
calculations  plotted  in  Figure  17  using  Eqs.  (2)  and  (13)  over  a  range  of  equivalent  ice 
sphere  radii  re  from  0.1  pm  to  350  pm.  This  comparison  shows  that  for  particle  sizes 
exceeding  ~  0.6  pm,  the  calculation  for  backscatter  ratio  (10.6/0.532)  using  scattering 
theory  for  spheres  lies  within  the  bounds  of  the  measured  backscatter  ratio  (10.6/0.532) 
from  cirrus  on  9  August. 
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Time  (GMT) 

Figure  13.  Contour  Intervals  of  -5.8,  -5.5,  -5.0,  and  -4.5  (m'1  sr‘‘)  on  a  Plot  of  log,0  (pn) 
at  0.532  nm  for  9  August  2001 


Time  (GMT) 


Figure  14.  Contour  Intervals  of -8.5,  -8.0,  -7.5, -7.0,  and  -6.5  on  a  Plot  of  log10(/?Jr)at  10.6 
|^m  for  9  August  2001 


The  same  technique  is  also  used  for  determining  yn  for  1 .064  pm  and  0.532  pm 
backscatter.  Corresponding  to  the  logio  ( P„)  contours  plotted  in  Figure  15,  the  measured 
backscatter  ratio  yn  (1.064/0.532)  is  plotted  in  Figure  18  as  a  function  of  time  and 
altitude.  The  measured  ( 1 .064/0.532)  backscatter  ratio  range  for  9  August  is 
approximately  0.5  to  1 .30.  The  theoretical  calculation  of  y n  as  a  function  of  equivalent 

radius  re  for  ice  spheres  is  plotted  in  Figure  19  along  with  the  range  of  measured  yn . 
Except  for  the  region  between  20  pm  and  100  pm,  the  measured  range  of  the  ratio  yn 

encloses  the  theoretical  values  for  equivalent  particle  sizes  exceeding  ~  2.0  pm. 

However,  for  either  case,  yn  (10.6/0.532)  or  yn  (1.064/0.532),  little  information  can  be 

gleaned  concerning  a  preferred  equivalent  or  mode  cirrus  particle  size  by  comparing  the 
measured  range  of  values  and  the  theoretical  calculations  over  the  range  of  equivalent  ice 
particle  radii. 

To  better  understand  the  representativeness  of  these  results,  the  measured  backscatter 
ratio  calculations  are  repeated  for  two  more  cirrus  measurement  days,  14  and  16  August. 
The  logio  P„  contour  analyses  for  14  August  are  shown  in  Figures  20  through  22  for  the 
three  wavelengths.  As  seen  from  these  figures  and  Figure  7,  the  cirrus  altitudes  on  these 
days  run  from  ~  9.5  -  12  km  with  the  principal  cirrus  layer  centered  near  10  km.  The 
range  of  log10  )  magnitudes  at  10.6  pm  was  comparable  to  those  measured  on  9 

August,  -  8.5  m'1  sr*1  to  -  6.5  m'1  sr*1,  but  cloud  center  values  of  log10  (/?„)  were  slightly 

lower  for  X  =  0.532  pm  and  X  =  1.06  pm,  with  maximum  values  of-  5.0  and  -  5.5 
observed,  respectively.  The  typical  measured  y ^  ratio  range  for  (10.6/0.532)  for  14 
August  was  1  x  10'3  to  5  x  10'2  according  to  Figure  23.  This  encompassed  a  lower 
particle  size  limit  for  the  theoretical  calculation  of  ~  re  =  0.5  pm  but  exceeds  by  a  factor 
of  5  the  maximum  value  of  calculated  yn  as  depicted  in  Figure  17.  On  14  August,  the 
range  of  measured  backscatter  ratio  yn  (1.064/0.532)  at  cirrus  altitudes  is  about  0.1  to  0.5, 
according  to  Figure  24.  This  encompasses  only  those  particles  whose  equivalent  radius  is 
less  than  ~  2.0  pm  according  to  Figure  19. 
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Backscatter  Ratio 


Figure  17.  Backscatter  Coefficient  Ratio  yn  for  Ice  Spheres  as  a  Function 
of  Equivalent  Radius  for  X=10.6  pm  ,  0.532  pm.  Dashed  lines  represent 
range  of  measured  yn  in  cirrus  for  9  August. 
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Figure  19.  Backscatter  Ratio  7^  for  Ice  Spheres  as  Function  of  Equivalent 
Radius  for  X=1.064  pm  ,  0.532  pm.  Dashed  lines  represent  range  of 
measured  y,  values  in  cirrus  for  9  August  2001 
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Figure  22.  Contour  Intervals  of  -6.5,  -6.0,  -5.5,  -5.0,  and  -4.5  (m'W)  on  a  Plot 
of logj0  ( )  at  X  =  1 .06  |nn  for  14  August 


Figure  23.  Contour  Intervals  of  0.003,  0.01,  0.1  and  0.3  on  a  Plot  of 
yK  (10.6/0.532)  for  14  August 
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Figure  24.  Contour  Intervals  of  0.01 , 0.25, 0.5,  and  1 .0  on  a  Plot  of 
y„  (1.064/0.532)  for  14  August  2002: 


Time  (GMT) 

Figure  25.  Contour  Intervals  of  -6.0,  -5.5,  -5.0,  and  -4.5 
on  a  Plot  oflogI0(/?;r)  at  k  =  0.532  |am  for  16  August  2001 
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Figure  26.  logl0  (y^)  at  X  =  10.6  fxm  for  16  August  2001:  Contour 
Interval:  -8.5, -8.0, -7.5, -7.0  (m'1  sr'1) 
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Figure  27.  log10  at  X  =  1 .064  jam  for  1 6  August  200 1 :  Contour 
Intervals:  -6.5,  -6.0,  -5.5,  -5.0,  -4.5 
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Figure  28.  ^(10.6/0.532)  Analysis  for  16  August  2001 :  Contour  Interval 
0.003,0.01,0.10 


Figure  29.  yn  (1 .064/0.532)  for  1 6  August  2001 :  Analysis  for  Contour  Intervals 
0.1.  0.5.  1.0.  1.5 


Figure  30.  Backscatter  Ratio  y„  for  Long  Perpendicularly-oriented  Ice  Cylinders 
a  Function  of  Cylinder  Radius  for  X=10.6  um ,  0.532  um 
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Figure  31 .  Backscatter  Ratio  yn  for  Long  Perpendicularly-oriented 
Cylinders  as  a  Function  of  Cylinder  Diameter  for  1=  1.064  urn  ,  0.532 


For  16  August,  the  measured  cirrus  altitudes  spanned  the  region  from  about  9  km  to  12.  5 
km,  encompassing  a  temperature  range  of-30C  to  -60C  (see  Figures  8  and  1 1).  The 
measured  log]0  (yS^)  time  vs.  height  contour  plots  at  cirrus  altitudes  for  each  of  the  three 

wavelengths  are  depicted  in  Figures  25  through  27.  log10  (/l^)  magnitudes  ranged  from 

-8.5  to  -  6.5  m’1  sr'1  for  X  =  10.6  pm  and  from  -  6.0  to  -  4.5  m'1  sr'1  for  0.532  pm  with 
cloud  center  magnitudes  exceeding  -  7.0  m’1  sr’1  and  -  4.5  m'1  sr’1  for  10.6  pm  and  0.532 
pm,  respectively.  The  range  of  measured  lidar  cirrus  backscatter  log10  (/?,,)  for  all  three 
days  is  summarized  by  wavelength  in  Table  4. 

Table  4.  Measured  Lidar  Cirrus  Backscatter  Coefficient  log10  (/?, )  Ranges  for  9, 14 
and  16  August  2001 

Date  0.532  um  1.06  urn  10.6  urn  Altitude  Range _ Temp.  Range 

9  Aug  -6.0  to -4.5  -6.0  to -5.0  -8.5  to -6.5  12.8  km- 15.3  km  -55Cto-70C 

14  Aug  -6.0  to -5.0  -6.5  to -5.5  -8.5  to -6.5  9.0  km -13.0  km  -30Cto-60C 

16  Aug  -6.0  to —4.5  -6.5  to —5.5  -8.5  to -6.5  9.0  km—  13.5  km  -30Cto— 60C 

The  measured  yK  (10.6/0.532)  in  cirrus  for  16  August  ranged  from  0.003  to  0.01  (Figure 

28),  encompassing  the  theoretical  ratio  curve  for  particle  sizes  exceeding  0.6  pm  (see 
Figure  17),  except  for  near  the  center  of  the  cirrus  layer  where  values  at  times  equal  or 
exceeded  0.10.  The  measured  ^(1.064/0.532)  ratio  varied  from  0.1  to  ~  1.0  before  18 
GMT  and  stayed  roughly  within  the  0.1  to  0.5  range  after  18  GMT  (Figure  29).  The 
yn (1.064/0.532)  measurements  before  18  GMT  roughly  enclose  the  ^ratios  expected 
from  theory  in  figure  19  but  the  measurements  after  18  GMT  are  generally  at  the  low  end 
of  the  theoretical  backscatter  ratio  values,  indicative  of  particle  sizes  less  than  a  couple  of 
pm  in  radius. 

The  range  of  yn  ratios  calculated  from  lidar  measurements  and  scattering  theory  for  ice 
spheres  (1  pm  <  re  <  100  pm)  are  summarized  in  Table  5. 
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Table  5.  Measured  Lidar  Cirrus  Backscatter  Coefficient  Ratio  Ranges  for 
9, 14  and  16  August  2001 


Date  y„  (10-6/0.5321 

9  Aug  0.003-0.01 
14  Aug  0.001-0.05 
16  Aug  0.003-0.10 
Theory  0.004-0.01 


y„  ( 1 .06/0.532)  Altitude  Ranee _ Temp.  Ranee 

0.5-1.30  12.8  km- 15.3  km  -55Cto-70C 

0. 1  -  0.50  9.0  km  -  1 3.0  km  -30C  to  -60C 

0.1  -1.00  9.0km- 13.5  km  -30Cto-60C 

0.3-2.00 


It  is  also  useful  to  compare  the  lidar-measured  cirrus  with  the  ^computed  from  theory 
for  long  ice  cylinders  or  rods  using  Eq.  (11).  The  theoretical  yn  ratios  (10.6/0.532)  and 
(1.064/0.532)  for  ice  cylinders  are  plotted  in  Figures  30  and  31,  respectively,  along  with 
the  range  of  measured  yn  for  the  9  August  cirrus  layers.  The  computed  yn  (10.6/0.532) 
for  ice  cylinders  exceeds  that  for  spheres  by  almost  an  order  of  magnitude  over  the  range 
of  equivalent  radii  while  the  computed  yn  (1.064/0.532)  for  ice  cylinders  exhibits 
magnitudes  similar  to  that  calculated  from  spheres.  Thus,  the  measured  range  of 
yn  (10.6/0.532)  for  9  August  only  accounts  for  a  relatively  small  size  range  of  equivalent 
radii  (<  1.0  pm)  theoretically  associated  with  long  ice  cylinders.  However,  the  measured 
range  of  ^(1.064/0.532)  encloses  nearly  the  entire  range  of  cylinder  radii  associated 
with  calculated  yn  from  0.1  pm  to  0.1  mm.  Thus  for  either  case,  using  spheres  or 
cylinders,  no  significant  information  can  be  deduced  concerning  preferred  particle  size. 
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This  last  conclusion  is  also  supported  by  the  results  contained  in  Appendix  C  which 
describe  the  normalized  backscatter  and  extinction  integration  kernel  (scattering 
efficiency  per  particle  radius)  behavior  as  a  function  of  ice  sphere  radius  for  the  three 
lidar  wavelengths. 
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7.  ESTIMATES  OF  ICE  WATER  CONTENT  FROM  CIRRUS  BACKSCATTER 

The  bulk  microphysical  properties  of  cirrus  can  be  characterized  by  particle  number 
density  N0  or  ice  water  content  piwc .  If  the  particle  size  distribution  n{r)  is  known  for  a 

polydispersion  of  ice  spheres,  the  ice  water  content  piwc  can  be  computed  as: 

Pi»c=^fpice[nir)r3dr  (kg  m'3)  (14) 

where  pice  is  the  density  of  ice  (which  is  a  weak  function  of  temperature).  Using  Eqs. 

(2),  (12),  and  (14),  the  measured  backscatter  and,  assuming  an  equivalent  or  modal 
sphere  radius,  piwc  can  be  computed. 

Analogously,  for  a  polydispersion  of  ice  cylinders,  the  ice  water  content  is  expressed  as 

Piwc  =  *  Pice  [n  (0 1  r2dr  (kg  m'3 )  (15) 

where  /  is  the  cylinder  length.  Note  that  the  inversion  of  Eq.  (1 1)  needed  for  equation 
(15),  using  Eq.  (12),  yields  number  density  N0  times  cylinder  length  / ,  N0-l  (m~2) , 

instead  of  number  density  N0  (m~3)  as  is  the  case  with  spheres  (Eq.  (2)). 

Calculations  of  time  series  of  piwc  for  a  polydispersion  of  ice  spheres  for  the  lidar 

measurements  on  9  August  at  14  km  altitude,  14  August  at  10  km  altitude,  and  16  August 
at  10  km  altitude  are  shown  in  Figures  32  through  34.  Results  are  displayed  assuming 
equivalent  radii  of  10  pm  and  100  pm.  The  particular  altitudes  selected  represent  the 
location  near  where  maximum  f5n  was  measured  in  the  lower  half  of  the  cirrus  altitude 

regions.  For  9  August,  using  the  measured  backscatter  at  14  km  altitude  and  temperature 
of  -65C,  the  time  series  of  calculated  ice  water  contents  show  relatively  good  agreement 
between  wavelengths  for  ice  water  content  magnitudes  exceeding  1  x  10'3  g  m'3 . 

Typical  ice  water  content  values  calculated  ranged  between  1  x  10'3  and  5  x  10'2  g  m'3  for 
re  =  100  pm  and  between  2  x  1  O'4  to  3  x  1 0'3  g  m'3  for  re  =  10  pm.  The  calculations  are 
compared  with  other  values  listed  in  the  literature.  Heymsfield  (1977)  reported  ice  water 
content  values  for  cirrus  at  a  temperature  of  -50C  ranging  from  ~  3  x  10"4  g  m  ‘3  to  8  x  10' 
3  g  m'3  using  particle  size  spectra  measurements  from  in-situ  particle  measurement 
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probes.  Heymsfield  (1994)  reported  a  mean  ice  water  content  of  3.9  x  1  O'3  g  m'3  for 
tropical  cirrus  in-situ  measurements  made  near  Kwajalein  with  temperatures  between  - 
60C  and  -70C.  This  value  is  plotted  as  a  dashed  line  on  Figure  32. 

Time  series  calculations  of  ice  water  content  are  plotted  in  Figures  33  and  34  using  the 
backscatter  measurements  at  10  km  altitude  on  14  August  and  16  August,  respectively, 
for  all  three  wavelengths.  On  both  days,  the  air  temperature  at  10  km  was  approximately 
-35C.  Heymsfield  ( 1 994)  listed  ice  water  contents  computed  from  in-situ  particle  size 
probe  measurements  in  cirrus  as  a  function  of  temperature  made  during  the  FIRE 
campaign  and  averaged  over  a  1  km  interval.  For  temperatures  near  -35C,  measured  ice 
water  contents  for  cirrus  not  associated  with  deep  synoptic-scale  systems  averaged  8.6  x 
10'3  g  m'3  with  a  median  value  of-  1  x  10'2  g  m‘3.  This  mean  value  is  plotted  as  a 
dashed  line  on  the  figures,  /^calculations  for  14  August  show  these  magnitudes 

regularly  being  achieved  only  for  the  measurements  at  10.6  pm  assuming  re  =  100 
pm,  and  are  not  achieved  using  f3n  measurements  at  1 .06  pm.  For  1 6  August,  the  mean 

in-situ  ice  water  content  magnitudes  are  achieved  routinely  for  measurements  at  10.6  pm 
and  more  frequently  for  /^measurements  at  0.532  pm,  again  assuming  an  equivalent 

radius  of  100  pm.  In  no  instance,  for  either  14  or  16  August,  are  plW  magnitudes  of  1  x 
10'  gm'  achieved  for  Px  measurements  at  1.06  pm.  It  should  be  noted  that  the 

backscatter  signal  at  0.532  pm  is  self-calibrated  against  the  Rayleigh  component.  While 
the  same  technique  of  calibration  is  also  applied  to  the  1 .06  pm  signal,  a  weaker  Rayleigh 
component  typically  introduces  a  higher  uncertainty. 

Calculations  of  ice  water  content  were  also  accomplished  for  the  scattering  efficiencies  of 
long  circular  cylinders  on  9  August  at  14  km  altitude  for  the  three  wavelengths  and  are 
plotted  in  Figure  35.  The  ratios  of  ice  water  content  for  spheres  ( iwc )  and  ice  water 
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IWC  (g  nr3)  IWC  (g  m 3)  IWC  (g  nr3) 


9  August  2001  z  =  14  km  0.532  nm 


9  August  2001  z  =  14  km  1 .064  (im 


Figure  32.  Ice  Water  Content  (iwc)  Calculations  for  9  Aug  Using  Backscatter 
Measurements  and  Assuming  a  Modified  Gamma  Polydispersion  of  Ice  Spheres 


IWC  (g  m-J) 


14  August  2001  z  *  10  km  0.532  nm 


14  August  2001  z  =  1 0  km  1 .064  \im 


14  August  2001  z  =  10  km  10.6  (im 


Figure  33.  Ice  Water  Content  (g  m'3)  calculations  for  14  Aug  Using  Backscatter 
Measurements  and  Assuming  a  Modified  Gamma  Polydispersion  of  Ice  Spheres 
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IWC  (g  nr3) 


1 6  August  2001  z  =  1 0  km  0.532  |im 


Hour  (GMT) 


16  August  2001  z  =  10  km  1.064 


Figure  34.  Ice  Water  Content  (g  m 3)  Calculations  for  16  Aug  Using  Backscatter 
Measurements  and  Assuming  a  Modified  Gamma  Polydispersion  of  Ice  Spheres 
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9  August  2001  z  =  14  km  1.064 


17.0  18.0  19.0  20.0  21.0  22.0 

Hour  (GMT) 

9  August  2001  z  =  14  km  10.6  pm 


17.0  18.0  19.0  20.0  21.0  22  0 
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Figure  35.  Ice  Water  Content  (iwc)  Calculations  Using  Backscatter  Measurements 
and  Assuming  a  Polydispersion  of  Long  Circular  Cylindrical  Scatterers 


IWC  (g  nr3) 


9  August  2001  z  =  14  km  0.532  pm 


Figure  36.  Time  Series  of  piwc  Computed  for  9  August  at  14  km  From  the 
Derived  Measurements  for  A,  =  0.532  pm  for  the  Two  Equivalent  Radii 
10  pm  and  100  pm. 


9  August  2001  z  =  14  km  0.532  pm 


Figure  37.  The  Ratio  of  Ice  Water  Content  Calculated  using  Extinction  and 
Backscatter  Measurements  for  9  Aug  at  A,  =  0.532  pm  for  the  Two  Equivalent 
Radii  10  pm  and  100  pm  is  Plotted  as  a  Function  of  Time. 


content  for  long  circular  cylinders  ( iwcr )  are  also  listed  by  the  figures  for  each  equivalent 

radius  re.  For  the  wavelengths  of  0.532  jam  and  1.064  pm,  the  ratio  is  less  than 

iwcr 

one,  ~  0.3  to  0.4.  However,  for  10.6  pm,  the  ratio  is  significantly  greater  than  one, 

iwcr 

being  ~  7.  This  implies,  that  for  a  given  particle  equivalent  radius,  the  integrated 
backscatter  efficiency  Qbscl  is  significantly  larger  for  long  circular  ice  cylinders  than  for 

ice  spheres  at  X.  =  10.6  pm.  For  X  =  0.532  pm  and  1.064  pm,  however,  integrated 
Qbsd  f°r  ice  spheres  is  2  -  3  times  the  size  of  the  backscatter  efficiency  for  ice  cylinders. 

Using  Eq.  (14)  and  Eq.  (19)  below,  the  ice  water  content  is  also  calculated  using  the 
derived  extinction  values  for  X  =  0.532  pm  that  are  presented  in  the  next  section.  A 
time  series  of  the  ice  water  content  values  calculated  using  the  extinction  coefficients  is 
plotted  for  9  August  using  the  two  equivalent  radii  of  10  pm  and  100  pm  in  Figure  36. 
The  ratio  of  ice  water  content  computed  from  backscatter  measurements  and  ice  water 

content  computed  from  the  extinction  measurements  VWCp*  is  also  plotted  for  the  two 

iwcp 

Pext 

equivalent  radii  and  displayed  in  Figure  37.  The  ice  water  content  ratio  varies  from  ~ 
0.55  for  re=  10  pm  to  ~  0.80  for  re  =  100  pm.  The  proximity  of  the  calculated  ratios  to 
unity  indicate  a  similarity  in  the  numerical  integration  for  backscatter  and  extinction  at  X 
=  0.532  pm  when  the  entire  particle  size  distribution  is  spanned.  This  similarity  exists 
despite  the  non-spherical  shape  of  the  measured  ice  particles  and  the  distinct  difference  in 
character  of  backscatter  and  extinction  efficiencies  Qbscl  and  Qat  at  optical  wavelengths  as 
a  function  of  size  parameter  x  (see  Appendix  C). 

8.  RANGE-RESOLVED  EXTINCTION  AT  LIDAR  WAVELENGTHS  IN 
CIRRUS 

The  extinction  coefficient  for  the  Nd:YAG  lidar  wavelength  0.532  pm  can  be  calculated 
using  a  technique  previously  applied  to  high  altitude  cirrus  by  Nee  et  al  (1998).  In  this 
technique,  it  is  assumed  that  the  extinction-to-backscatter  ratio 
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S  where 


s=^44  h  <i6> 

PAZ) 

is  an  empirical  constant  within  the  cloud  and  p^,  is  the  wavelength-dependent  extinction 
coefficient.  For  a  given  cloud,  an  effective  value  of  S  can  be  determined  using  the 
equation 

s=2'n(ffe))[C/?'<z)^"1  lsr)  <17) 

where  Z\  and  Z2  are  the  lower  and  upper  altitude  limits  of  the  cloud.  The  value  f5„  (z2)  is 


essentially  the  backscatter  measured  at  the  top  of  the  cloud  and  PK  (z,)  is  the  backscatter 
at  the  base.  Since  backscatter  outside  the  cloud  is  dominated  by  Rayleigh  scattering  at 
cirrus  altitudes,  the  upper  value  pn  (z2)  is  set  to  the  value  computed  from  Rayleigh 
scattering  theory.  The  Rayleigh  volume  backscatter  coefficient  for  air  is  expressed  as 


9  m  - 1 
NX*  m2  +2 


(m-V'1) 


where  m  is  the  refractive  index  of  air,  N  is  the  number  of  air  molecules  per  unit  volume, 


and  X  is  the  optical  wavelength  (0.532  pm) .  Here  N  =  — —  where  N A  is  Avogadro’s 

M 

number,  p  is  the  density  of  air,  and  M  is  the  molecular  weight  of  air.  A  derivation  of 
Eq.  (18)  is  presented  in  Appendix  D. 


Figure  38  shows  a  vertical  profile  of  the  Rayleigh  volume  backscatter  coefficient  Pn  for 
the  two  lidar  wavelengths,  X  —  0.532  pm  and  X  =  1.06  pm  computed  from  the 
radiosonde  launch  on  9  August  (see  Figure  9).  For  0.532  pm,  the  characteristic  Pn  value 

is  6  x  10'7  m'1  sr'1  at  10  km  and  3  x  10"7  m"1  sr'1  at  15  km.  For  X=  1.064  pm,  the 
characteristic  PK  value  is  4  x  10'8  m’1  sr'1  at  10  km  and  2  x  10'8  m'1  sr'1  at  15  km. 


For  9  August  and  16  August  data,  the  extinction-to-backscatter  ratio  S  was  determined  to 
be  50  and  28,  respectively,  from  the  integrated  lidar  backscatter  measurements  at  X- 


0.532  pm.  For  14  August,  a  value  of  S  =  30  was  used,  consistent  with  the  ratio  used  by 
Nee  et  al  (1998)  and  consistent  with  the  observations  by  Sassen  and  Cho  (1992).  The 
extinction-to-backscatter  ratio  S  can  also  be  calculated  as  a  function  of  the  particle  size 
using  scattering  theory  along  with  the  expression  for  (Eq.  (2))  and  the  expression  for 
extinction  by  single  scattering  by  a  polydispersion  n(r) : 

Pex,  =  f  Qex, n  {r)nr2dr  (nT1)  (19) 

Qext  *s  ^e  extinction  efficiency  for  ice  spheres,  which  is  a  function  of  the  wavelength, 
particle  size,  and  particle  index  of  refraction.  Here,  the  extinction  efficiency  for  spheres 
is  defined  as  a  function  of  the  scattering  amplitude  Sl  evaluated  at  scattering  angle  0  =  0: 


where  5,  (0)  is  evaluated  after  van  de  Hulst  (1957)  as 

5>(0)  =  7Z(2"  +  1)k+O  (21) 

Z  n= 1 

The  extinction  efficiency  for  uniform  perpendicularly-oriented  long  cylinders  is 
determined  in  an  analogous  form  to  that  described  by  Eq.  (8)  for  cylindrical  backscatter. 
Results  of  the  theoretical  calculations  for  S  for  both  spheres  and  ice  cylinders  are 
displayed  in  Figure  39  for  a  range  of  ice  sphere  equivalent  radii  from  0.1  pm  to  350  pm 
along  with  the  measured  range  of  values  (28  -  50).  For  particle  radii  greater  than  ~1  pm, 
the  theoretical  calculations  of  S  for  ice  spheres  range  between  18  and  28,  slightly  less 
than  that  observed  from  the  lidar  backscatter  measurements.  Values  of  S  calculated  for 
ice  cylinders  with  equivalent  radii  greater  than  ~  10  pm  lie  within  the  range  of  lidar- 
measured  values. 

The  altitude-dependent  extinction  coefficient  is  calculated  as 

Pe*,(z)  =  S-PAz)  (m'')  (22) 

where  S  is  assumed  to  be  constant  within  the  cloud.  There  are  provisions  to  recursively 
correct  the  backscatter  to  eliminate  the  effect  of  extinction.  For  these  profiles,  the 
corrected  backscatter  is  then  used  to  recalculate  the  extinction  profile.  For  most  cases, 
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the  optical  thickness  is  sufficiently  small  so  that  no  recursive  correction  is  applied.  In 
general,  the  error  associated  with  the  evaluation  of  Eq.  (22)  is~lxl05m1  (0.01  km  '). 
There  are  two  systematic  biases  associated  with  the  derived  extinction  profile.  First,  S  is 
applied  to  all  altitudes  over  the  cloud  layer  depth  including  those  with  and  without  cloud. 
This  may  result  in  an  overestimate  of  extinction  since  S  is  much  smaller  (~  8ti/3)  for 
Rayleigh  scatterers  than  for  the  cloudy  part.  Another  source  of  bias  is  due  to  the  solar 
background  subtraction  and  “shutter  function”  effect.  These  types  of  bias  should  not  be  a 
problem  if  the  higher  error  level  of  10'5  m'1  is  accepted  for  the  extinction. 

Time  vs.  altitude  contour  plots  of  extinction  coefficient  (3ext ,  computed  for  A,  =  0.532  pm 

using  the  measured  backscatter  and  equations  described  above,  for  9, 14,  and  16  August, 
are  shown  in  figures  40  through  42,  respectively.  Near  the  center  of  the  cirrus  layer  from 
13.8  km  to  14.5  km  around  20  GMT  on  9  August,  the  extinction  coefficient  at  0.532  pm 
exceeds  0.5  km'1 .  Over  a  2  km  layer  depth,  the  extinction  coefficient  exceeds  0.1  km'1 
for  a  period  on  the  order  of  one  hour  on  9  August.  On  14  August,  the  cirrus  extinction 
coefficient  is  less  than  on  9  August  with  values  exceeding  0.1  km'1  only  near  the  center  of 
the  cirrus  layer  ~  10  km.  For  16  August,  some  layers  approaching  a  2  km  depth  over 
which  the  extinction  coefficient  exceeds  0.5  km  1  are  observed  around  1830  GMT  near  10 
km  and  some  areas  exceeding  0.3  km 1  are  found  at  higher  altitudes  between  12  and  13 
km  after  20  GMT.  These  results  indicate  that  higher  altitude  cirrus  (~  14  km) 
unconnected  with  any  distinct  weather  system  or  thunderstorm  activity  can  be  associated 
with  as  large  or  larger  optical  extinction  values  (>  0.5  km  *)  than  those  associated  with 
cirrus  at  lower  altitudes  (~  9  km). 

For  comparison  with  passive  remote  sensing  retrievals,  the  cirrus  optical  depth  tc  is 
calculated  using  the  altitude-dependent  lidar-measured  extinction  coefficient  for  X  = 
0.532  pm.  This  is  plotted  as  a  time  series  in  Figure  43  for  the  three  measurement  days. 
Here,  the  cirrus  optical  depth  rc  is  defined  as 

=  t  Pe*{z)dz  (23> 

JZ, 
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where  Zx  and  Z2  represent  the  altitude  integration  limits  for  the  observed  cirrus  layers. 
The  dashed  lines  in  Figure  43  denote  the  optical  depth  regions  associated  with  subvisual 
cirrus  (rc  <  0.03),  thin  cirrus  (0.03  <  tc  <  0.30),  and  cirrus  up  to  the  attenuation-limited 


Rayleigh  Backscatter  from  Air  Molecules 


Figure  38.  Rayleigh  Volume  Backscatter  Coefficient  (5n  for  Air  Molecules  for  9 
August  2001  19  GMT  Hanscom  AFB 
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S  (SI) 


532  nm  Extin ction/Backscatter  Ratio 


r„(nm) 


Figure  39.  Extinction  to  Backscatter  Ratio  S:  Measured  vs.  Theoretical 
Calculations  of  S  at  X  =  0.532  pm  for  a  Polydispersion  of  Ice  Spheres  (Dashed 
Line)  and  Long  Ice  Cylinders  with  Light  at  Perpendicular  Incidence  (Solid 
Line!. 
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0.532  pm  Cirrus  Optical  Depth 


17.0  18.0  19.0  20.0  21.0  22.0  23.0 

Time  (UTC) 


Figure  43.  Time  Series  of  Measured  Cirrus  Optical  Depth  tc(X  =  0.532  pm)  for  9 

August,  14  August,  and  16  August  2002.  Numbers  in  parentheses  denote  maximum 
observed  values.  Dashed  lines  mark  descriptive  regions  of  cirrus  type  according  to 
optical  depth  (subvisual,  thin,  and  semi-opaque)  after  Sassen  and  Cho  (1992). 
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opacity  (0.3  3.0)  (Sassen  and  Cho,  1992).  However,  the  range  for  subvisual 

values  lies  below  the  optical  depth  error  threshold  here  for  a  10  km  thick  layer  (~  0.10). 

For  9  August,  the  integration  limits  Z,  and  Z2  are  10  km  and  16  km,  respectively,  while 
for  14  and  16  August,  the  integration  limits  are  6.6  km  and  17  km.  For  9  August,  the 
mean  cirrus  optical  depth  rc  is  0.39  with  a  maximum  value  of  1 .40  and  minimum  value  of 

0.15.  For  14  August,  the  mean  cirrus  optical  depth  rcis  ~  0.35  with  a  maximum  of  0.56 

and  minimum  of  0.06.  For  16  August,  the  mean  cirrus  optical  depth  is  ~  0.29  with  a 
maximum  value  of  1 .86  and  minimum  value  of  0.02.  Excluding  the  peak  areas  in  Figure 
43,  and  using  the  optical  depth  criteria  of  Sassen  and  Cho  (1992),  thin  cirrus  is  the 
predominant  type  during  the  9  August  and  16  August  measurements  (prior  to  20  GMT), 
while  semi-opaque  cirrus  is  the  characteristic  cirrus  type  for  14  August. 

9.  SUMMARY  AND  CONCLUSION 

Range-resolved  lidar  backscatter  measurements  of  cirrus  were  made  at  three  wavelengths 
-  0.532  pm,  1.064  pm,  and  10.591  pm  -  over  a  temperature  range  of-30C  to  -70C  and 
altitude  range  of  9  km  to  15.5  km  on  several  days  in  August  2001  at  Hanscom  AFB, 
Massachusetts.  Three  days  were  selected  for  analysis  and  study  -  9  August,  14  August, 
and  16  August.  The  cirrus  backscatter  measurements  were  analyzed  in  combination  with 
scattering  theory  calculations  to  better  understand  their  behavior  with  respect  to 
wavelength,  characteristic  particle  size,  ice  water  content,  and  comparative  optical 
properties. 

Lidar-measured  volume  backscatter  coefficients  fiK  were  analyzed  for  the  three 
wavelengths  as  a  function  of  time  and  altitude.  For  X  =0.532  pm  and  1.064  pm,  logt0 
Pn  magnitudes  ranged  from  -  4.5  to  -  6.0  m'1  sr'1  within  the  cirrus  cloud  area,  with 
values  of -5.0  frequently  occurring  near  the  cloud  center,  while  log)0  Pn  magnitudes  for 
X.  =  10.6  pm  ranged  typically  from  -  6.5  to  -  8.5  m'1  sr1  with  typical  cloud  center  values 
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near  -7.0  .  The  measured  volume  backscatter  coefficient  ratios  yK  were  determined  as  a 
function  of  time  and  altitude  for  the  wavelength  ratios  10.6/1.064  and  10.6/0.532. 

The  measured  backscatter  ratio  magnitudes  yn  were  compared  with  magnitudes  computed 
from  scattering  theory  for  both  ice  spheres  and  long  circular  cylinders  using  a  modified 
gamma  particle  size  distribution  «(r) .  Measured  ^(10.6/0.532)  during  the  three 

measurement  days  ranged  from  1  x  10’3  to  1  x  10'1  while  measured  ya( 1.064/0.532) 

ranged  from  0.1  -  1.30.  General  agreement  existed  between  the  overall  range  of 
measured  backscatter  ratio  magnitudes  and  those  computed  from  theory  for  both 
wavelength  pairs,  but  there  was  no  consistent  information  which  could  be  discerned 
concerning  characteristic  cirrus  particle  size  associated  with  the  ^measurements  and  the 
behavior  of  the  theoretical  backscatter  ratio  function  as  a  particle  size. 

Assuming  the  modified  gamma  particle  size  distribution  and  given  an  equivalent  or 
modal  particle  size,  calculations  of  the  time  series  of  ice  water  content  piwc  were  made 
from  the  backscatter  measurements  at  the  three  wavelengths  using  scattering  theory. 
Typical  values  of  piwc  computed  for  9  August  20  GMT  at  14  km  ranged  from  ~1  x  10 3  g 

m'3  to  ~5  x  10‘2  g  m'3  for  an  equivalent  sphere  radius  of  100  pm  while  the  calculated  ice 
water  content  for  an  equivalent  sphere  radius  of  10  pm  was  smaller  by  more  than  an 
order  of  magnitude.  />im.  values  significantly  exceeding  1  x  10'2  g  m'3 ,  given  re  =  100 
pm,  were  observed  on  16  August  for  for  both  X  =  0.532  pm  and  10.6  pm  and  for  10.6 

pm  on  14  August.  Similar  calculations  were  also  performed  for  long  ice  cylinders 
assuming  uniform  orientation  and  light  of  perpendicular  incidence.  For  the  shorter 
wavelengths  (0.532  pm  and  1.064  pm),  piwc  calculated  assuming  cylinders  exceeded  that 

for  spheres  by  a  factor  of  2  or  3,  while  piwc  calculated  assuming  spheres  exceeded  that 
computed  for  ice  cylinders  by  a  factor  of  seven  for  backscatter  at  10.6  pm. 


Using  knowledge  of  the  Rayleigh  backscatter  coefficient  (m'1  sr'1)  for  air  molecules,  the 
cirrus  backscatter,  and  assuming  a  constant  value  in-cloud,  estimates  of  extinction-to- 
backscatter  ratio  S  (sr)  in  cirrus  were  computed  for  9  and  16  August  for  X  =  0.532  pm. 
The  measured  values  of  S  ranged  from  28  sr  to  50  sr.  The  measured  values  of  S  agreed 
generally  with  theoretical  values  for  both  ice  spheres  and  long  ice  cylinders  for  particle 
size  distributions  where  the  equivalent  radii  exceeded  ~  1  pm.  The  derived  values  of  S 
were  applied  to  the  time  and  height  profiles  of  the  measured  cirrus  to  produce 

corresponding  estimates  of  cirrus  particle  extinction  (3ext  (m"1).  These  estimates  showed 
/L  at  X  =  0.532  pm  exceeding  0.5  km'1  on  9  August  and  16  August  near  14  km  and  near 
10  km  altitude,  respectively,  and  exceeding  0.1  km'1  in  the  vicinity  of  10  km  on  14 
August.  These  cirrus  extinction  profiles  correspond  to  peak  optical  depth  values  zc  of 
1.4  and  1.9  on  9  August  and  16  August,  respectively,  and  0.6  on  14  August. 

These  results  also  imply  that  accurate  determination  of  characteristic  and  bulk 
microphysical  properties  of  cirrus  clouds  will  be  difficult,  at  best,  using  only  lidar 
backscatter  measurements.  This  is  a  result  of  the  nearly  constant  backscatter  coefficient 
ratio  y  K  for  optical  wavelength  pairs  over  a  wide  range  of  modal  or  characteristic  particle 
size  parameters.  For  cirrus,  the  characteristic  sizes  of  the  polydispersion  usually  exceed 
the  few  pm  threshold  where  the  onset  of  a  nearly  constant  lidar  backscatter  ratio  yn  is 

encountered.  The  availability  of  backscatter  measurements  at  longer  wavelengths, 
such  as  X  =  3.2  mm  or  X  =  8.6  mm  for  cloud  radar,  should  significantly  increase  the 
variation  of  ys  with  size  parameter  when  paired  with  optical  wavelengths. 
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APPENDIX  A.  CALCULATION  OF  SCATTERING  COEFFICIENTS  FOR  A 
SPHERE 


A  rather  straightforward  algorithm  for  calculating  the  spherical  scattering  coefficients  and 
handling  their  recursion  problem  was  derived  by  Aden  (1951).  Here  the  scattering 
coefficients  are  defined  as 


Jn(*)  fg.M-Mg.OO] 

h}2){x)\pn{x)-man(y)\ 


Jjx)  \<rm(y)-m<Tm(x)} 

K[2){x)\(Tn{y)-mpn(x)\ 


where 


and 


A.-i(2)(*)  n 

h}2)(x)  x 


Ws^ln[^W]= 


J'n- 1  (g) 
Jn{x) 


n 

x 


are  expressed  in  terms  of  the  Riccati-Bessel  functions  jn  (x),nn  (x) . 


(Al) 


(A2) 


(A3) 


(A4) 


For  the  lowest  orders,  integers  n  =  0,  n=  1 
xjn  (x) ,  xnn  (x)  are  defined  as: 

xj0  (x)  =  sinx 

.  ,  x  sinx 
XJ\  (*)  = - cosx 


(A5) 
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xw0(x)  =  -cosx 


XM, 


-cosx 

- sinx 

X 


(A6) 


Noting  that  y  -  mkr  -mx  ;  m  =  m[A)  =  mr  -  imi ,  the  complex  index  of  refraction  where 
i  =  7-T  ,  requires  the  use  of  complex  variables  such  that 

z  =  x  +  iy  then  \zf  =  (x2  +/)  ;  also  z=x-iy,  \z\  =  \z\ 

Given  y  =  mrx-imix  then 

sin  y  =  sin  (mrx)  cosh  (-m,x) + i  cos(wrx)sinh  (-/w,x) 
cos  y  =  cos  (mrx)  cosh  (-m,x)  -  i  sin  ( mrx )  sinh  (-m(x) 

Higher  order  terms  of  the  Riccati-Bessel  functions  can  be  determined  using  recursion 
relations  for  spherical  Bessel  functions  —also  known  as  Bessel  functions  of  fractional 
order: 


A*.  (x)  =  ^y^7.(x)-y_,(x) 


(A7) 


and  h}2)(x)  =  jn(x)-inn(x) 

which  is  the  second  Hankel  function  of  order  n.  The  Hankel  function  employs  the  same 
recursion  relations  as  other  spherical  Bessel  functions. 

In  the  case  of  spherical  Bessel  functions,  if  significant  imaginary  components  exist  in  the 
index  of  refraction,  the  arguments  for  the  hyperbolic  trigonometric  functions  can,  as  x 
increases,  become  too  large  for  the  machine  data  type  and  floating  point  overflow  will 
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occur.  This  can  be  avoided  by  expressing  (4)  in  exponential  form  and  factoring  out 
exponential  terms  in  the  numerator  and  denominator  which  have  a  positive  exponent. 

The  recursion  relation  for  the  logarithmic  derivatives  derived  from  Infeld  (1947)  is: 


x2  +  nxan_A  (x)  -  n2 
nx-x2crn_i  (x) 


(A8) 


Thus,  using  the  logarithmic  derivatives  (3,4),  recursion  relations,  and  (8),  the  scattering 
coefficients  a„  and  b„  can  readily  be  calculated. 

Given  the  expression  for  scattering  efficiency  Qsct  and  extinction  efficiency  Qexl  (van  de 
Hulst,  1957) 


X  n= 1 


(A9) 


2e*=4Z(2M+1)Rek+6'') 

X  n= 1 


(A10) 


convergence  criteria  for  evaluating  all  sphere  scattering  efficiency  calculations  were  set 
arbitrarily  as: 


|G«,  (")-&*("-  0|y  ^  | G«fW-G«r("-1)| 


xr_ 

2 


<10" 


Stephens  (1961)  discussed  the  numerical  mode  which  develops  when  calculating  the 
extinction  and  scattering  efficiencies  Qexl ,  Qsct  for  water  spheres  and  becomes  negatively 
unbounded  as  the  size  parameter  x  — >  0.  Encounters  with  this  mode  may  occur,  in 
practice,  at  microwave  or  lower  frequencies.  One  remedy  is  to  use  van  de  Hulst’s  (1957) 
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series  expansion  for  small  x  (~  x  <  0.01)  in  lieu  of  the  standard  determination  using 
amplitude  coefficients: 


Qex,=~  Im 


Q,„=x4 


4x 


8 

m2-\ 

3 

m2  +2 

m  -1  4  3 

— x - +— x3 

m2  +  2  15 

2 

+  ... 


^m2-lV 
m2  +  2 


mA  +  27  m2  +38 
2m2  +  3 


>  +  x4  Re  < 

8 

2  ' 

> 

l3 

^m2  +2, 

>  +  ... 


(All) 


For  very  small  x  and  mx  Rayleigh  scattering  (e.g.  x  <  0.01),  series  expansions  should  also 
be  used  for  calculation  of  the  scattering  amplitudes  S}  and  S2.  Following  van  de  Hulst 
(1957),  these  are: 


isx  r 

=  —^~ [3  +  (3t  +  3wcos#  +  5wcos#)x2  -3isx3  +  ...J  (A12) 

isx^  r 

S2  =  —  [3cos 6  +  (3t cos 9  +  3«  +  5wcos  26) x2  -  3 isx3  cos 6  + ...]  (A1 3) 


where 


2  m2  -1  _3  m2  —2 

3  m2 +2’  1  5m2  +  2’ 


1  m2  +2 
w  = - - - ; 

10  w2+3 
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APPENDIX  B:  CALCULATION  OF  SCATTERING  EFFICIENCIES  FOR  A 
LONG  CIRCULAR  CYLINDER  (ROD) 

Solutions  for  the  scattering  coefficients  an  and  b„  at  perpendicular  incidence  are 
expressed  in  terms  of  the  Bessel  functions  of  integer  order  (van  de  Hulst  (1957),  Kerker 
(1969)): 


b  _  (>Q j \ M K_ (*) 

mJ„'  (y)Hn(x)-Jn(y)Hn'  (x) 


(Bl) 


Jl(y)Jn(x)-mjn(y)j’(x) 
Jn  ( y)  Hn  (a:)  -  mJ„  (y)  H'n  (x) 


(B2) 


Here,  J„  is  the  Bessel  function  of  the  first  kind  of  order  n  where  primes  denote  first 
derivatives,  m  is  the  complex  index  of  refraction  m  =  mr  -imi  of  the  scatterer  where 

i  =  4-\ ,  x  = - ,  the  size  parameter,  and  y  =  mx.  H„(x)  is  the  Hankel  function  of  the 

A 

second  kind,  Hf*  (x)  =  J„  ( x)  -iY„  (x) ,  where  Y„(x)  is  the  Bessel  function  of  the  second 
kind  of  order  n.  Y„(x)  is  also  referred  to  as  the  Neumann  or  Weber  function. 

The  Bessel  functions  of  integer  order  are  calculated  using  a  series  expansion  for  x  <  30 
and  an  asymptotic  expression  for  x  >  30.  For  Bessel  functions  of  the  first  kind  Jn,  these 
expressions  are: 

Jv(z)  =  ±{~l)  r —  ( series )  x<30  (B3a) 

vK  ’  n\  T(v  +  «  +  l) 


JAz)= 

Z~*  00 


nz  cos 


z — V7i — n 


( asymptotic )  x  >  30 


(B3b) 
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where  v  and  n  are  integers,  T  is  the  gamma  function,  and  z  may  be  real  (x)  or  complex 

oo. 


For  Bessel  functions  of  the  second  kind  Y„,  the  Weber  or  Neumann  functions,  the 
equations  used  are: 


+  l)+y'('1 +  *  +  !))  7:7-— L:  ( series )  x<30 


(B4a) 


x  Ho v  '  '  '  '  "  £!(«  +  £)! 


^(l)  =  ~y  (Euler's  constant  ~  0.5772....) 


where 


^(«)  =  -r+Xr' 


y„  (z)  =  yj2/(xz) sin (z-\nn -\n)  {asymptotic)  x>30  (B4b) 

2-*co 

Equations  (B3)  and  (B4)  can  be  found  in  Chapter  9  of  Abramowitz  and  Stegun  (1964) 
and  other  texts  such  as  Arfken  (1985).  Appropriate  series  and  asymptotic  size  parameter 
regimes  were  determined  empirically  and  checked  against  the  Harvard  computation 
tables  (1947),  tables  in  Abramowitz  and  Stegun  (1964),  and  Jahnke  and  Emde  (1945). 

In  practice,  only  the  functions  for  integer  orders  0  and  1  are  determined  using  the  given 
values  of  x  andy.  The  higher  order  Bessel  functions  needed  for  evaluation  of  Eqs.  (Bl) 
and  (B2)  including  the  first  derivatives  are  calculated  using  the  following  recursion 
relations: 


forward: 

n.,(z)=yi:(*)- n-,6) 

(B5a) 

backward: 

-,„-i  (z)=^-J.(z)-J.Az) 

(B5b) 
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first  derivative: 


J'n{z)  =  J„- 1(2) — -4(z) 

Z 


(B5c) 


Forward  recursion  can  be  used  for  Bessel  functions  of  the  second  kind  (Weber  functions) 
while  backward  recursion  must  be  used  for  Bessel  functions  of  the  first  kind.  The 
backward  or  downward  recursion  is  performed  using  a  technique  described  by 
Abramowitz  and  Stegun  (1964)  and  attributed  to TCP.  Miller  (Bickley  et  al,  1952).  As 
opposed  to  the  spherical  or  Bessel  functions  of  fractional  order,  the  appearance  of 
numerical  error  in  forward  recursion  using  Eq.  (B5a)  appears  for  integer  orders  as  small 
as  9  or  10.  Also,  the  downward  recursion  algorithm  seems  to  require  more  care  in 
determining  the  “starting”  n  value  compared  to  the  spherical  Bessel  functions. 


A  hybrid  recursion  scheme  was  used  due  to  the  floating  point  overflow  encountered  in 
the  asymptotic  form  of  Eq.  (B3b),  which  occurs  when  a  significant  imaginary  component 
mi  *  x  is  reached  for  the  complex  size  parameter  y.  Using  “standard”  procedures,  the 
scattering  coefficients  are  expressed  explicitly  in  terms  of  the  Bessel  functions.  Another 
approach  is  to  rewrite  the  scattering  coefficients  in  terms  of  the  derivative  of  the 
logarithm  for  the  Bessel  functions  of  the  first  kind  with  complex  argument,  and  by 
factoring,  eliminate  the  positive  imaginary  exponents  that  cause  the  overflow.  To 
accomplish  this,  Eqs.  (Bl)  and  (B2)  are  re-written  as 


b  _  Mg,  00 /,(*)- 4,' (*) 

m(Tn{y)Hn(x)-Hn'(x) 


(B6) 


a„  = 


<r*{y)Jn(x)~mJ«  (x) 

n(y)HAX)~mHn'{X ) 


(B7) 


where 


Infeld’s  (1947)  recursion  relation 


°.(y)  = 


ny-y2<r„- 1  (y) 


(B8) 
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is  applied  to  the  logarithmic  basis  functions  (orders  0  and  1)  and  using  downward 
recursion,  the  scattering  coefficients  can  then  be  calculated  for  any  order  n.  Given 
y  =  xr+  xi  where  xr  =  mrx  and  x(.  =  -mtx ,  the  derivative  of  the  logarithmic  basis 
function  for  order  0,  using  Eq.  (B3b),  is: 


a  /  v  l  +  [cos(2xr)-fsin(2xf)](-/g2^) 
°  ^  -/+  cos^x^-isin^x,)]^ 


(B9) 


Evaluation  of  Eqs.  (B3)  -  (B9)  were  accomplished  for  size  parameters  ranging  from  1  x 
10'10  to  100,000  and  for  integer  orders  exceeding  80,000.  The  Bessel  function 
calculations  were  compared  with  tables  9.1  and  9.4  in  Abramowitz  and  Stegun  (1964)  for 
size  parameters  x  and  integer  orders  up  to  100  and  with  the  tables  of  Chapter  8  in  Jahnke 
and  Emde  (1945). 


Given  the  expression  for  scattering  efficiency  Qsct  and  extinction  efficiency  Qat  (van  de 
Hulst,  1957) 


Case  1 :  v  =0 ,  (E  ||  axis) 


Rebo  +  2XRe6» 

n=\ 


n= 1 


(B10) 


Case  2:  u  =  0 ,  (H  ||  axis) 


o  oo 

&,,=-£  Rea, 


2 

x 


Ren0  + 


2lR=a 


+2£kf 


n=\ 


(B 11) 
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where  E  and  H  denote  the  electric  and  magnetic  field  vector  amplitudes,  respectively, 
convergence  criteria  for  evaluating  all  cylindrical  scattering  efficiency  calculations  were 
set  arbitrarily,  using  the  Case  1  functions,  as: 


and  |0ett(»)-att(«-l)|- 


<  10- 
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APPENDIX  C:  NORMALIZED  BACKSCATTER  AND  EXTINCTION  KERNEL 
FUNCTIONS  FOR  ICE  SPHERES  AT  LIDAR  WAVELENGTHS 

Given  expressions  for  volume  backscatter  and  extinction  coefficients  for  a  polydispersion 
of  spheres: 

, 6k  ^bsct  n  (r)Kr2dr  (m'1  sr'1)  (Cl) 

Pm  =  f  Qex,  (m>  r>  £)  n  (r) nrl<ir  (m4)  (C2) 


these  can  be  written  in  a  generalized  form 

g(x,A,m)  =  pK(x,A,m)f(x)dx  (C3) 

*Xa 


where  jc  s  lnr ,  g  (x,A,m)  is  either  the  volume  backscatter  or  extinction 
coefficient,  K (x,  2,  m)  is  a  kernel  function,  and  /  (x)  is  the  size  distribution  function.  In 
practice  (e.g.  Qing  et  al,  1989),  f(x)  can  be  expressed  as  the  non-dimensional  volume 
density  function 


f(x) 


4  n  3  dN 
Tr  dlnr 


(C4) 


given  n(r)  =  ,  and  K  is  expressed  as: 


K(x,A,m)=  —  Ql  (m'1  sr'1  or  m'1  ) 
4  r 


(C5) 


where  QL  = 


or  Qexr 


Normalized  plots  of  the  kernel  functions  K  as  a  function  of  ice  sphere  radius  for  the  three 
lidar  wavelengths  are  shown  below. 


For  lidar  backscatter  kernels  associated  with  ice  spheres,  the  particle  radius  range  of 
influence  runs  from  ~  0.1  pm  to  ~  10  pm  for  the  shorter  wavelengths  (0.532  pm  and 
1.064  pm)  and  ~  1  pm  to  ~10  pm  for  X  =  10.6  pm.  For  extinction  kernels,  the  shorter 
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lidar  wavelengths  are  characterized  by  a  relatively  narrow  but  smoother  radius  range  of 
influence  (0.1  (am  -  5.0  (am)  compared  to  backscatter.  Extinction  at  the  10.6  pm 
wavelength  retains  a  strong  influence  for  particle  sizes  up  to  10  pm  and  this  influence 
does  not  decrease  until  the  sphere  radius  exceeds  10  pm  and  diminishes  in  a  manner 
similar  to  the  other  wavelengths.  The  markedly  different  behavior  at  X  =  10.6  pm, 
compared  to  the  shorter  wavelengths,  is  likely  due  to  the  large  imaginary  component  of 
the  refractive  index  for  ice  and  its  strong  component  of  ice  particle  absorption. 


Backscatter  Kernels  for  Ice  Spheres 


Figure  Cl .  Kn  (r,  A, m)  (normalized) 

for  wavelengths  X  =  0.532  pm,  1.064 
pm,  and  10.6  pm. 


Extinction  Kernels  for  Ice  Spheres 


Figure  C2.  Kal(r,A,m)  (normalized) 

for  wavelengths  X  =  0.532  pm,  1.064 
pm,  and  10.6  pm. 
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APPENDIX  D:  RAYLEIGH  BACKSCATTER  CROSS-SECTION  AND  VOLUME 
BACKSCATTER  COEFFICIENT  FOR  AIR 

From  van  de  Hulst  (1957),  the  backscatter  cross-section  or  is 

cr  =  47tr2  l(r,0  =  7c)/lo  ( m 2)  (Dl) 

where  I0  is  the  intensity  of  the  incident  light,  /  is  the  intensity  of  the  scattered  light  at  a 
large  distance  r  from  the  scattering  particle  and  9  is  the  scattering  angle. 

For  isotropically  polarized  Rayleigh  scattering,  the  scattering  intensity  is 

(l  +  cos2#)&4|a|2/0 

/=i - £ r 44-  <W> 

Thus 

<r  =  2;r(l  +  cos20)&4|  a  |2  where  (D3) 


k  =  —  and  a  is  the  polarization  volume,  a  is  expressed  using  the  Lorentz-Lorenz 

A 

formula  applicable  to  molecular  optics  (van  de  Hulst  (1957)): 


1  3(m2-l) 

4 jtN  [m2  +2^ 


(m3) 


(D4) 


Here  N  =  ^~ 
M 


the  number  of  molecules  per  irnit  volume,  where  NA  is  Avogadro’s 


number,  p  is  the  density  of  air,  and  Mis  the  molecular  weight  of  air.  In  Eq.  (D4),  m  is  the 
refractive  index  of  air,  which  is  a  real  quantity.  Substituting  Eq.  (D4)  into  Eq.  (D3) 
yields 


cr  =  • 


8  tzN'- 


-(l  +  cos2#)A:4 


m2  - 1 


m2  +2 


(D5) 
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After  Liou  (1980),  the  average  scattering  cross-section  ex  is  expressed  in  terms  of  the 
scattering  coefficient  J3scl  (m1)  and  N  as 

a  =  Psc,/N  (D6) 

If  the  scattering  angle  6  =  n,  the  volume  backscatter  coefficient  can  be  written  as 


9k4 

m2- 1 

2_  9;r2 

m1- 1 

4  K 

16;r2iV 

w2  +  2 

~  NX4 

w2  +2 
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